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CHRISTIAN MERCAT 

Abstract. Wc detail the theory of Discrete Ricmaim Surfaces. It 
takes place on a cellular decomposition of a surface, together with 
its Poincare dual, equipped with a discrete conformal structure. 
A lot of theorems of the continuous theory follow through to the 
discrete case, we will define the discrete analogs of period matrices, 
Riemann's bilinear relations, exponential of constant argument and 
series. We present the notion of criticality and its relationship with 
integrability. 
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1. Introduction 

Riemann surfaces theory was a major achievment of XlXth century 
mathematics, setting the framework where modern complex analysis 
bloomed. Nowadays, surfaces are intensively used in computer science 
for numeric computations, ranging from visualization to pattern recog- 
nition and approximation of partial differential equations. A lot of 
these computations involve, at the continuous level, analytic functions. 
But very few algorithms care about this specificity, although analytic 
functions form a relatively small vector space among the space of func- 
tions, problems are usually crudely discretized in a way that doesn't 
take advantage of analyticity. 

The theory of discrete Riemann surfaces aims at filling this gap and 
setting the theoretical framework in which the notion of discrete ana- 
lyticity is set on solid grounds. 

Most of the results in this paper are a straightforward application 
of the continuous theory [H [2] together with the results in [31 HI [5], 
to which we refer for details. We define the discrete period matrix, 
which is twice as large as in the continuous case: the periods of a 
holomorphic form on the graph and on its dual are in general different, 
but the continuous limit theorem, given a refining sequence of critical 
maps, ensures that they converge to the same value. The main tool is 
the same as in the continuous case, the Riemann bilinear relations. 

2. Discrete Riemann surfaces 

2.1. Discrete Hodge theory. We recall in this section basic defini- 
tions and results from [1] where the notion of discrete Riemann surfaces 
was defined. We are interested in discrete surfaces given by a cellular 
decomposition <0> of dimension two, where all faces are quadrilaterals 
(a quad-graph [61 [7J [8]). Its vertices and diagonals define, up to ho- 
motopy and away from the boundary, two dual cellular decompositions 
T and T*: The edges in r* are dual to edges in Ti, faces in T% are 
dual to vertices in T and vice-versa. Their union is denoted the double 
A = T U r*. A discrete conformal structure on A is a real positive 
function p on the unoriented edges satisfying p(e*) = l/p(e). It defines 
a genuine Riemann surface structure on the discrete surface: Choose 
a length S and realize each quadrilateral by a lozenge whose diagonals 
have a length ratio given by p. Gluing them together provides a flat 
riemannian metric with conic singularities at the vertices, hence a con- 
formal structure [9]. It leads to a straightforward discrete version of 
the Cauchy-Riemann equation. A function on the vertices is discrete 
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Figure 1 . The vertices and diagonals of a quadrilateral 
define a pair of dual edges. 




ke* 



0. The vertex dual to a face. 



1. Dual edges. 




2. The face dual to a vertex. 



Figure 2. Duality. 

holomorphic iff for every quadrilateral (x,y,x',y') G ^>2) 
(2.1) f(y') - f(y) = ip(x,x') (f(x') - f(x)) . 



We recall elements of de-Rham co homology, doubled in our context: 
The complex of chains C(A) = Co (A) © Ci(A) © C 2 (A) is the vec- 
tor space span by vertices, edges and faces. It is equipped with a 
boundary operator d : Cfc(A) — > Cfe_i(A), null on vertices and fulfilling 
d 2 = 0. The kernel ker d =: Z,(A) of the boundary operator are the 
closed chains or cycles. Its image are the exact chains. It provides the 
dual spaces of forms, called cochains, C k (A) := Hom(Cfc(A), C) with a 
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Figure 3. The discrete Cauchy-Riemann equation. 
coboundary d : C k (A) — > C fc+1 (A) defined by Stokes formula: 



df := / (d(x, x')) = f(x') - f(x), J J da:= jot. 

(x,x>) F dF 

A cocycle is a closed cochain and we note a e Z k (A). 

These spaces are equipped with the canonical scalar product, weigthed 
according to p on edges and averaged on the graph and its dual: 

(a,/?):=±X>(e)(7a > ) ( 1(5 
Duality of complexes allows us to define a Hodge operator * on forms 

by 

* : C k (A) -> C 2 " fe (A) 
C°(A)9/ ^ */ : JJ *f.-f(F*), 

(2.2) C^A) 9 a i-> *a : /*<*:= -p(e*) / a, 

*/ e J e* 

C 2 (A) 9w h-> *u; : (*a>)(x) := / / u. 



It fulfills * 2 = (— ld C k) k . The endomorphism A := — d * d* — * d* d 
is the usual discrete Laplacian: Its formula on a function at a vertex 
x G T with neighbours x±, . . . , xy £ T is the usual weighted averaged 
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difference: 

v 

(A(/))(x) = ^p(x,x fc )(/(x)-/(x fc )). 
k=i 

The space of harmonic forms is defined as its kernel . 

The Hodge star and the laplacian are real operators. Since * 2 = —Id 
on functions, it is natural to consider them on complexified cochains. 
The discrete holomorphic forms are special complex harmonic forms: 
a 1-form 

(2.3) a G C 1 (A) is holomorphic iff da = and * a = —ia, 

that is to say if it is closed and of type (1,0). Let a", resp. d" the 
compositions of the exterior derivative with the projection on the space 
of (1,0), resp. (0, l)-forms, Eq. (I2.3P is equivalent to d'a = 0. We will 
note a G fi 1 (A). A function / : Ao — » C is holomorphic iff df is 
holomorphic, which is equivalent to (12. ip and we note / G fi°(A). 

In the compact case, — * d* is the adjoint d* of the coboundary 
operator d and the Hodge theorem orthogonally decomposes forms into 
exact, coexact and harmonic, 

C k (A) = Imd© 1 Im d* ® L Ker A, 

harmonic forms are the closed and co-closed ones, and harmonic 1-form 
are the orthogonal sum of holomorphic and anti-holomorphic ones: 

Ker A = Ker d D Ker d* = Ker d' 0" 1 Ker d" . 



2.2. Wedge product. We construct a wedge product on <0> such that 

• the canonical weighted hermitian scalar product reads as ex- 
pected 

(a,p) = JJ aA *p, 

• and the coboundary operator d$ on <), is a derivation for this 
product A : C7 fc «» x C7^) -> C k+l (<». 
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It is defined by the following formulae, for /, g £ C°(^)), a,(3 E C 1 ({)) 
and £ C 2 (<»: 

(2.4) (/■ <7)(z) :=/(*) -<7(z) forx£^ , 

/(*) + f(y) 



(2.5) 



/ • a := 



a for £ Oi, 



(2.6) 

(xi,X2,X3,Xi) 

(2.7) ///.. 

(xi,X2,X3,X4) 



^ : =iE / a J p- j a 

fc_1 (x fc _ 1 ,x fe ) (xj,,x fc+ i) (x fc+1 ,x fc ) (x fe ,x fe _i) 
/(xi)+/(x a )+/(x 3 )+/(ai4) 



(xi,X 2 ,X 3 ,X4) 

for (xi,X2,x 3 ,x 4 ) £ <> 2 - 



A form on (} can be averaged into a form on A: This map A from 
C*(0) to C*(A) is the identity for functions and defined by the following 
formulae for 1 and 2-forms: 



/ 



(2.8) 



(2-9) 



(x,x') 



+ + 



\(x,y) (y,x') {x,y') (y',x') J 



fc=l 



(x k ,y k ,x,y k _ 1 ) 



where notations are made clear in Fig. HI The map A is neither injective 
nor surjective in the non simply-connected case, so we can neither define 
a Hodge star on <0> nor a wedge product on A. Its kernel on 1-forms is 
Ker (A) = Vect (d$e), where e is the biconstant, +1 on V and —1 on Y*. 
But d^A = A d^ so it carries cocycles on <0 to cocycles on A. Its image 
are these cocycles of A verifying that their holonomies along cycles of 
A only depend on their homology on the combinatorial surface. Given 
a 1-cocycle \x £ Z 1 (A) with such a property, a corresponding 1-cocycle 
v £ 2' 1 (<0) is built in the following way: Choose an edge (£0,2/0) £ Oij 
for an edge (x, y) £ <0>i with x and x on the same leaf of A, choose 
two paths X X)XQ and X yo>y on the double graph A, from x to x and y 
to y respectively, and define 



(2.10) 



v :-- 



VQ'V 



[-] 
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where [7] = [A^,, + (xo, Vo) + \ ,y + (y, %)] is the class of the full cycle 
in the homology of the surface. Changing the base points change fi by 
a multiple of d$e. 

It follows in the compact case that the dimensions of the harmonic 
forms on <) (the kernel of AA) modulo de, as well as the harmonic 
forms on A with same holonomies on the graph and on its dual, are 
twice the genus of the surface, as expected. Unfortunately, the space 
Im A = TL 1 - © Im d is not stable by the Hodge star *. We could never- 
theless define holomorphic 1-forms on <) but their dimension would be 
much smaller than in the continuous, namely the genus of the surface. 
Criticality provides conditions which ensure that the space *Im A is 
"close" to Im A. 



y 

.0. 



X 



X 



xy' 

y 




flUD 



to 



Figure 4. Notations. 



We construct an heterogeneous wedge product for 1-forms: with 
a, /3 e C7 1 (A), define a A (3 e C^O) by 



(2.H) 



a A (3 :- 



a J J J 

\(x,x') (y,y>) (y,y>) (x',x) ) 



\ 



It verifies A(a$) A A((3$) = ot^ A /?<>, the first wedge product being 
between 1-forms on A and the second between forms on (}. The usual 
scalar product on compactly supported forms on A reads as expected: 



(2.12) w) = §£p(e)QQ Ui) = J J 

eGAl 6 6 02 



a A *(3 
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2.3. Energies. The L 2 norm of the 1-form df, called the Dirichlet 
energy of the function /, is the mean of the usual Dirichlet energies on 
each independant graph: 

(2.13) 

E D {f):=\\\df\\ 2 = \{df,df) = \ Y, P(x,x')\f(x')-f(x)\ 2 

(x,x')eA! 

E D (f\r) + E D (f\ r .) 
2 

Harmonic maps minimize this energy among functions fulfilling certain 
boundary conditions. 

The conformal energy of a map measures its conformality defect, it 
is null on holomorphic functions: 

(2.14) E c U):=\\\df-i*df\\ 2 . 

It is related to the Dirichlet energy through the same formula as in the 
continuous: 

E c {f) = \{df-i*df,df-i*df) 

= l\\df\\ 2 + \\H*df\\ 2 + ±Re(df, -i*df) 

= l\\df\\ 2 + flm ff dfAdf 



2 



(2.15) =E D (f)-A(f) 

where the area of the image of the application / in the complex plane 
has the same formula 



(2.16) A(f) = % // dfAdf 

as in the continuous case since, for a face (x, y, x', y') G O2, the algebraic 
area of the oriented quadrilateral ( fix), f(x'), f(y), f(y') ) is given by 



df A df = i Im ({f(x') - f(x))(f(y') - f(y)) 
= -2iA(f(x),f(x'),f(y),f(y')). 



3. Period matrix 



We use the convention of Farkas and Kra [T] , chapter III, to which we 
refer for details. Consider (<),p) a discrete compact Riemann surface. 
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3.1. Intersection number, on A and on For a given simple (real) 
cycle C G Zi(A), we construct a harmonic 1-form rjc such that § A T]c 
counts the algebraic number of times A contains an edge dual to an 
edge of C: It is the solution of a Neumann problem on the surface cut 
open along C (see [3] for details). It follows from standard homology 
technique that r\c depends only on the homology class of C (all the 
cycles which differ from C by an exact cycle OA) and can be extended 
linearly to all cycles as 77. : -Hi (A) — > C 1 (A); it fulfills, for a closed form 

e, 

(3.1) <fd= [[ n c A6, 

Jc J Jo 

and a basis of the homology provides a dual basis of harmonic forms 
on A. Beware that if the cycle C G Zi(T) is purely on T, then this 
form 7]c\ r = is null on T. 

The intersection number between two cycles A,Be Z±(A) is defined 

as 

(3.2) A - B :— JJ^VaAVb. 

It is obviously linear and antisymmetric, it is an integer number for 
integer cycles. Let's stress again that the intersection of a cycle on T 
with another cycle on T is always null. A cycle C G Zi((}) defines 
a pair of cycles on each graph Cr G Zi(F), Cr* G Zi(T*) which are 
homologous to C on the surface, composed of portions of the boundary 
of the faces on A dual to the vertices of C. They are uniquely defined 
if we require that they lie "to the left" of C as shown in Figj5l By 
the procedure (I2.10p applied to r]c r + r]c r , , we construct a 1-cocycle 
Tj C G Z 1 ^) unique up to de, and since WO, de A 9 = 0, Eq. (13. 2p defines 
an intersection number on Zi(§). Unlike the intersection number on A, 
this one has all the usual expected properties. In particular Eq. (13.21) 
holds for A,Be Zi(0). 

3.2. Canonical dissection, fundamental polygon. The complex 
<0 being connected, consider a maximal tree T C (}\, that is to say T 
is a Z 2 -homologically trivial chain and every edge added to T forms a 
cycle. A canonical dissection or cut-system K of the genus g discrete 
Riemann surface <0> is given by a set of oriented edges (ek)\<k<2 g such 
that the cycles K C (TUet) form a basis of the homology group Hi((}) 
verifying, for 1 < k, i < g 



(3.3) N fc -N* = 0, K fc+9 -K £+9 = 0, N fc -N /+fl = £ M . 
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Figure 5. A path C on <0> defines a pair of paths Cr 
and Cr* on its left. 

They actually form a basis of the fundamental group 7Ti(<0>) and the 
defining relation among them is (noted multiplicatively) 

(3-4) n» fc Nrf^ = l. 

k=l 

The construction of such a basis is standard and we won't repeat the 
procedure. What is less standard is the interpretation of Eq. (13 .4|) in 
terms of the boundary of a fundamental domain, discretization intro- 
duces some subtleties (that can safely be skept in first instance). We 
end up with the familiar 2g x 2g intersection numbers matrix on ^>. 

Considering T U 6fc clS 8b rooted graph, we can prune it of all its 
pending branches, leaving a simple closed loop K^, attached to the 
origin O by a simple path X k (see Fig. E]), yielding the cycle These 
three cycles are deformation retract of one another, ^ C ^ C T U 
hence are equal in homology. 

In the continuous basis of the homology can be realized by 

2g simple arcs, transverse to one another and meeting only at the base 
point. It defines an isometric model of the surface as a fundamental 
domain homeomorphic to a disc and bordered by Ag arcs to identify 
pairwise. In the discrete case, by definition, the set <0\N of the cellular 
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Figure 6. 1. A maximal rooted tree in a quadrilateral 
decomposition of the torus. 2. An additional edge defines 
a rooted cycle K 1; pruned of its dangling trees. 3. Its un- 
rooted version, the simple loop H^. 



complex minus the edges taking part into the cycles basis is homeomor- 
phic to a disc hence the surface is realized as a polygonal fundamental 
domain M. whose boundary edges are identified pairwise. 

But it is sometimes impossible to choose a basis of the homology 
verifying (I3.3P by simple discrete cycles which are transverse to one 
another. For instance, if the path is not empty, the cycle is 
not even simple. Moreover, some edges may belong to several cycles. 
In this case, the edges on the boundary of this fundamental polygon 
can not be assigned a unique element of the basis or its inverse, and 
therefore can not be grouped into only 4g continuous paths to identify 
pairwise but more than Ag. 

In fact, the information contained into the basis K is more than 
simply this polygon, the set of edges composing the concatenated cycle 

(3.5) Vi, K\ N^, N 2 , . . . , NJ 1 , K^ 1 ) 

encodes a cellular complex A4 + which is not a combinatorial surface 
and consists of the fundamental polygon M. plus some dangling trees, 
corresponding to the edges which belong to more than one cycle or 
participate more than once in a cycle (the paths A&), as exemplified 
in FigJTl By construction, the edge belongs to the cycle only, 
hence these trees are in fact without branches, simple paths whose 
only leaf is the base point O. To retrieve the surface, the edges of this 
structure AA+ are identified group-wise, an edge participating k times 
in cycles will have [k/2] +2 representatives to identify together, two on 
the fundamental polygon and the rest as edges of dangling trees. 

Eliminating repetition, that is to say looking at (13.51) not clS db SG~ 
quence of edges but as a simplified cycle (or a simplified word in edges), 
thins Ai + into Ad, pruning away the dangling paths. The fundamental 
polygon boundary loses its structure as 4g arcs to be identified pairwise, 
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Figure 7. Three different fundamental polygons of a 
decomposition of the torus (g = 1) by three quadrilat- 
erals: 1. The standard fundamental domain where the 
Ag paths are not adapted to {>. 2. .M+ is composed of 
edges of <0> composing 4g arcs (which may have portions 
in common) to identify pairwise, each edge corresponds 
to an element of the basis N or its inverse, except for 
edges of "dangling trees" which are associated with two 
such elements. 3. M. is composed of edges of <) compos- 
ing more than Ag arcs to identify pairwise, there is no 
correspondence with a basis of cycles. 



in general a basis cycle will be disconnected around the fundamental 
domain and a given edge can not be assigned to a particular cycle. 
This peculiarity gives a more complex yet well defined meaning to the 
contour integral formula for a 1-form 9 defined on the boundary edges 
of M+, 

(3.6) / 9 = V / 9 + I 9. 

This basis gives rise to cycles K r and N r * whose homology classes form 
a basis of the group for each respective graph, that we compose into 
N A defined by 

(3-7) N£ = C = 

^k+2g — ^fc+g' ^k+3g ~ n fc+g' 

for 1 < k < g so that while the intersection numbers matrix on {> is 
given by the 2g x 2g matrix 



(3.8) ■ N 



e)k,e 



/ 
-/ 
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the intersection numbers matrix on A is the 4g x 4g matrix with the 
same structure 



(3.9) (Nj • N 



r 






r 




(0 





I 


o\ 


r 











I 


p* 


-I 











p* 


{° 


-I 







r 



3.3. Bilinear relations. 

Proposition 3.1. Given a canonical dissection N, for two closed forms 

9,9>ez\0), 

(3.10) If 6 A6' = VM & I 6' - I 8 I 6' 

JJo j=1 i« J+9 Jn j+g Jkj t 

for two closed forms 6,6' G Z 1 (A), 

(3.11) [[ 6 A6' = y^ ( <[ 8 <[ 6' - I 6 I 6' 

J Jo UVk M + i s M +2B J «t , 

Proof [3A[ Each side is bilinear and depends only on the cohomol- 
ogy classes of the forms. Decompose the forms onto the cohomology 
basis (ojfe). On A, use Eq (13 . 1 5|) for the LHS and the duality property 
Eq. (13.141) for the RHS. On <), use their counterparts. (} 

Notice that for a harmonic form 6 G 7i 1 (A), the form *6 is closed as 
well, therefore its norm is given by 

(3.12) 6 G H\A) \\6\\ 2 = V ( & 6 I *6-<f 6 I *6 

3.4. Basis of harmonic forms, basis of holomorphic forms. We 

define a A , the basis of real harmonic 1-forms, dual to the homology 
basis K A , as described in Sec. 13. 1[ 

a k := %£ 10 an d 

v 

l k+2g 



k + 2g 

(3.13) a£ +2g := -t] h a for 1 < k < 2g 



which verify 



cue = ^ 



(3.14) & a l+2g 



kl, 



>kl, 



"fc + 2g 



14 CHRISTIAN MERCAT 

and dually, the intersection matrix elements are given by 



(3.15) 



K ■ K 



On <0>, the elements af 



defined up to de, verify A(af) 



V* k+3 and a% +g 



-rj^ k for 1 < k < g, 



a 



k+g. 



A(<*t+ 9 ) 



a 



and form a basis of the cohomology on <0> dual to K as well, 



k+2g 



a 



k+3g 







(3.16) 



a 



k+g 



k + g 



and 

for 1 < k < g, 



they fulfill the first identity in Eq. fl3.15p but the second is meaningless 
in general since * can not be defined on <). We will drop the mention 
A when no confusion is possible. 



Proposition 3.2. The matrix of inner products on A, 
(3.17) 



[Oik, Cttjkl 



1 < k < 2g, 
is a real symmetric positive definite matrix. 



A D 
B C 



Proof UTR It is real because the forms are real, and symmet- 
ric because the scalar product (12.121) is skew symmetric. Definition 
Eq. (I3.13P and Eq. (13.11) lead to the integral formulae. Positivity fol- 
lows from the bilinear relation Eq. (13.111) : for 9 = ^t=i£fc a ki with 
& G C, £? =1 |a| 2 > 0, 



2<? 

E 

0=1 

E 

k,t=l 
1.9 



*9 



2<? 



9 *9 

K 2 g+j J^2j 



3=1 



otk j *ae - 



Oik 



■29+3 



2j 



(3.18) 



^2 66(«fc,«f) > o. 



kl=l 



The form a k is supported by only one of the two graphs V or T*, 
the form *a k is supported by the other one, and the wedge product 
#p A 9' r = is null for two 1-forms supported by the same graph. 
Therefore the matrices A and C are g x g-block diagonal and B is 
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anti-diagonal. 
(3.19) 

'A T 



.4 



A r , 



D 





Bt.t* 



Bt*,t 




C 



C7 r , 








The matrices of intersection numbers ( 13. 9ft and of inner products differ 
only by the Hodge star *. Because * preserves harmonic forms and the 
inner product, we get its matrix representation in the basis a, 



—D 

-C 



A 
B 



B 



(3.20) 

and because * 2 
(3.21) 
(3.22) 
(3.23) 

On <0, while the Hodge star * can not be defined, we can obviously 
consider the following positive scalar product on the classes of closed 
forms modulo de, to which the set (a^) belong: 



C-A + I 
A-B 
C- l B 





t B ■ A 
B-C. 



[a 



(Aia^ 



t / f, ,a\ 
I J (x ,y) 



{x,y,x' ,y')e<>2 

p=p(x,x'),p*=p(y,y') 



I, 



(V ,x') 



J(x',y') 



a 



a 



\hv,x) a ) 



( ' +P+P* 


+p-p* 


-p- 


-P* 


-P+P*^ 


+p-p* 


+p+p* 


-PH 


-P* 


-p-p* 


-p-p* 


-p+p* 


+PH 


-P* 


+ P-P* 


y-p+p* 


—p—p* 


+P- 


-P* 


+P+P* J 



and it yields 

(3.24) (aJ,a?) W 



i?rr* + Bt*t 



>B 



'B T 



tt* t -or*r 



which, in general, can not be understood as the periods of a set of 
forms on <0> along the basis H. 

Let's decompose the space of harmonic forms into two orthogonal 
supplements, 

(3.25) n^A) = n\ ®^n\ 

where the first vector space are the harmonic forms whose holonomies 
on one graph are equal to their holonomies on the dual, that is to say 

(3.26) H\ := Vect (a k + a k+g , \<k<goi2g<k< 3g). 
Definition (13 . 1 3[) and Eq. (13. 1 P imply that 

(3.27) H\ = Vect (*a k - *a k+g , l<k<goi2g<k< 3g). 



ff(x, y) A 

I(y ,x') P 
f(x>,y>) & 
\I(y',x ) P) 



(3.28) 
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These elements in the basis (a k + ak+ g , ; ct k — Oi k+g ) for 1 < k < g and 
2g < k < 3g, are represented by the following invertible matrix: 

// 'Err* - *# r *r A r - A r » \ 

/ Cr — C*r* i?rr* — i?r*r 

'Err* + *5r*r A r + ^r* 

\0 C T + C r * Brr* + 5r*r/ 

It implies in particular that the lower right g x g block is invertible, 
therefore so is Eq. f)3.24p . 

3.5. Period matrix. 

Proposition 3.3. The matrix IT = C^ 1 ■ (i — B) is the period matrix 
of the basis of holomorphic forms 

(3-29) Ck:=(i-*)J2 C ki a ^ 
in the canonical dissection that is to say 

8 k ,i fori <k< 2g, 

Rk-2g,l for 2g < k < 4g, 



(3.30) <b Q 



k 



and IT is symmetric, with a positive definite imaginary part. 

The proof is essentially the same as in the continuous case pQ and 
we include it for completeness. 

Proof lcTR Let ujj := aj + i * ctj for 1 < j < Ag. These holomorphic 
forms fulfill 

(3.31) Pfcj := \{u k , Uj) = (a k , a,) + i (a k , - * atj) 

-i (' lOu, 1 < j < 2o 

(3.32) J *^° -J - y 



i ^k, 2g < j < Ag. 



P is the period matrix of the forms (u) in the homology basis K. The 
first 2g forms (^j)i<j<2g are a basis of holomorphic forms. It has the 
right dimension and they are linearly independent: 

2g 2g I 2g 2g \ 

+ ifJ>j){aj + i * Qfj) = ( ( A J + ^ k B ^) a i + Yl ^ k C i> k a2 9+i J 

j=l j=l \ k=l k=l J 

23 / 2g 2g \ 

(3.33) +i Y ( (a*j +Y^ k B 3J°) a i + Xk C i' k a2 9 + i ) 

j=l \ k=l k=l J 

is null, for A,/i G R only when A = /i = because C is positive 
definite. Similarly for the last 2g forms. The change of basis iC~ l on 
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them provides the basis of holomorphic forms (£). The last 2g rows of 
P is the 2g x Ag matrix (B — i I, C) hence the periods of (£) in K are 
given by (J, n). 

The first identity in Eq. fl3.30p uniquely defines the basis £ and a 
holomorphic 1-form is completely determined by whether its periods 
on the first 2g cycles of K, or their real parts on the whole set. 

Notice that because C is g x g block diagonal and B is anti-diagonal, 
II is decomposed into four g x g blocks, the two diagonal matrices form 
iC^ 1 and are pure imaginary, the other two form — ■ B and are 
real. 

Therefore the holomorphic forms (k are real on one graph and pure 
imaginary on its dual, 

(3.35) 1 < k < g => ( k e C^(T) © i C&(r*) 
g<k<2g (kEC^n^iC^T). 

We will call 

(3.36) n r = n r + Ik* 

the period matrix on the graph F the sum of the real periods of Ck, 
1 < k < g, on T, with the associated pure imaginary periods on the 
dual r*, and similarly for g < k < 2g, the period matrix on T*. 

It is natural to ask how close Hp and Hr. are from one another, and 
whether their mean can be given an interpretation. Criticality [31 H] 
answers partially the issue: 

Theorem 3.1. In the genus one critical case, the period matrices Hp 
and Ilr* are equal to the period matrix lis of the underlying surface S. 
For higher genus, given a refining sequence (() k ,Pk) of critical maps 
of E 7 the discrete period matrices Ilpfc and Il r *k converge to the period 
matrix lis • 

Proof WIS. The genus one case is postponed to Sec. 13.61 The con- 
tinuous limit comes from techniques in 0, H], developed in [S] which 
prove that, given a refining sequence of critical maps, any holomorphic 
function can be approximated by a sequence of discrete holomorphic 
functions. Taking the real parts, this implies as well that any harmonic 
function can be approximated by discrete harmonic functions. In par- 
ticular, the discrete solutions /& to a Dirichlet or Neumann problem on 
a simply connected set converge to the continuous solution / because 
the latter can be approximated by discrete harmonic functions and 
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the difference f k — gu being harmonic and small on the boundary, con- 
verge to zero. In particular, each form in the basis (ctf), provides a 
solution to the Neumann problem Eq. (13.161) and a similar procedure, 
detailed afterwards, define a converging sequence of forms (f, yielding 
the result. <0 

We can try to replicate the work done on A on the graph (}. A 
problem is that A r + A r and C-p + C-p need not be positive definite. 
Moreover, the Hodge star * doesn't preserve the space (A(a%)) of har- 
monic forms with equal holonomies on the graph and on its dual, so 
we can not define the analogue of a + i * a on (}. We first investigate 
what happens when we can partially define these analogues: 

Assume that for 2g < k < 3g, the holonomies of *a k on T are equal 
to the holonomies of *a k+g on T*, that is to say Cp = Cp* =: \C$ and 
Dpp* = Dr*r ='■ \D$. It implies that the transposes fulfill Byt* = 
Bp*p ='■ \B$ as well. We can then define f3 k _ g G Z l {(}) such that 
A(Pk-g) = * a k+gi uniquely up to de. The last g columns *(_B^,C^) of 
the matrix of scalar product Eq. f)3.24p are related to their periods in 
the homology basis in a way similar to Eq. (13.171) . By the same 
reasoning as before, the forms 



(3-37) C? = E C tk i « 9 ~ . l<k< 9 



=i 



verify A((j?) = ^ fc+g and have periods on ^ given by the identity 
for the first g cycles and the following g x g matrix, mean of the period 
matrices on the graph and on its dual: 

(3.38) U<> = C^-B,) = ^±^. 

The same reasoning applies when the periods of the forms *ak on 
the graph and on its dual are not equal but close to one another. In the 
context of refining sequences, we said that the basis (ocf), converges to 
the continuous basis of harmonic forms defined by the same Neumann 
problem Eq. (13.161) . Therefore 

(3.39) Cr-Cr* =o(l), B rr * - B r , r = o(l) . 

A harmonic form u k+g = o(l) on T* can be added to *a k + g such that 
there exists /3%_ g G Z l (<>) with A(/3^ g ) = *a k+g + v k+g , yielding forms 
Cjf , verifying A((£) = |(£ fe + ( k + g ) + o(l) and whose period matrix is 
11^ + o(l). Since the periods of a k converge to the same periods as its 
continuous limit, this period matrix converges to the period matrix lis 
of the surface. Which is the claim of Th. 13.11 
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In the paper [10], R. Costa-Santos and B. McCoy define a period 
matrix on a special cellular decomposition T of a surface by squares. 
They don't consider the dual graph T*. Their period matrix is equal to 
one of the two diagonal blocks of the double period matrix we construct 
in this case. They don't have to consider the off-diagonal blocks because 
the problem is so symmetric that their period matrix is pure imaginary. 



3.6. Genus one case. Criticality solves partially the problem of hav- 
ing two different g x g period matrices instead of one since they converge 
to one another in a refining sequence. However, on a genus one criti- 
cal torus, the situation is simpler: The overall curvature is null and a 
critical map is everywhere flat. Therefore the cellular decomposition is 
the quotient of a periodic cellular decomposition of the plane by two 
independant periods. They can be normalized to (1,t). The continu- 
ous period matrix is the 1 x 1-matrix r. A basis of the two dimensional 
holomorphic 1-forms is given by the real and imaginary parts of dZ on 
T and T* respectively, and the reverse. The discrete period matrix is 

the 2x2 matrix ( Jj 11 T T ) and the period matrices on the graph 



Re r Im r 

and on its dual are both equal to the continuous one. 

For illustration purposes, the whole construction, of a basis of har- 
monic forms, then projected onto a basis of holomorphic forms, yielding 
the period matrix, can be checked explicitely on the critical maps of the 
genus 1 torus decomposed by square or triangular /hexagonal lattices: 

Consider the critical square (rectangular) lattice decomposition of a 
torus <0> = (Tje 19 + 'Ze~ ie )/(2pe ie + 2qe~ ie ), with horizontal parameter 
p = tanO and vertical parameter its inverse. Its modulus is r = ^e 2lB . 
The two dual graphs T and T* are isomorphic. An explicit harmonic 
form a\ is given by the constant l/2p on horizontal and downwards 
edges of the graph V and on all the other edges. Its holonomies are 
1 and on the p, resp. q cycles. Considering l/2g and the dual graph, 
we construct in the same fashion a\, af*, a\* . The matrix of inner 
products is 



(3.40) 



{0i k ,0it) 



k.e 



1 

sin 26 



cos 



20\ 



cos 29 



cos 29 



\cos2# 
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using p+ 2 = l/sin26 l and p 2 = — 1/ tan 20 so that the period 
matrix is 

/o a-,\ tt 1 fi sin 20 cos 20 \ 

3.41 n = - _ . . _ . 

v ' p \ cos 20 i sin 20 y 

Therefore there exists a holomorphic form which has the same periods 
on the graph and on its dual, it is the average of the two half forms 
of Eq. f)3.30p and its periods are (1, ^e 2%e ) along the p, resp. q cycles, 
yielding the continuous modulus. This holomorphic form is simply the 
normalized fundamental form _ i6 . 

In the critical triangular /hexagonal lattice, we just point out to the 
necessary check by concentrating on a tile of the torus, composed of two 
triangles, pointing up and down respectively. We show that there exists 
an explicit holomorphic form which has the same shift on the graph 
and on its dual, along this tile. Let p~,p\ and p/ the three parameters 
around a given triangle. Criticality occurs when p_ p\ + p\ p/ + pj p_ = 
1. The form which is 1 on the rightwards and South- West edges and 
elsewhere is harmonic on the triangular lattice. Its pure imaginary 
companion on the dual hexagonal lattice exhibits a shift by % p\ in the 
horizontal direction and % (p\ + p_) in the North-East direction along 
the tile. Dually, on the hexagonal lattice, the form which is p\ p_ along 
the North-East and downwards edges and 1 — p\ p_ along the South- 
East edges, is a harmonic form. Its shift in the horizontal direction is 
1, in the North-East direction 0, and its pure imaginary companion on 
the triangular lattice exhibits a shift by i p\ in the horizontal direction 
and % (p\ + p_) in the North-East direction along the tile as before. 
Hence their sum is a holomorphic form with equal holonomies on the 
triangular and hexagonal graphs and the period matrix it computes is 
the same as the continuous one. This simply amounts to pointing out 
that the fundamental form dz can be explicitely expressed in terms of 
the discrete conformal data. 



4. Criticality and integrable system 

This theory can be viewed as the simplest (it is linear) of a series of 
integrable theories [Sj. We will present its quadratic counterpart, which 
leads to another version of discrete analytic functions, based on circle 
patterns. Along the way, we will see how discrete exponentials and 
discrete polynomials emerge due to integrable systems theory pieces of 
technology named the Backlund or Darboux transform [8]. 
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4.1. Criticality. Until now, everything has been purely combinato- 
rial, there was no reference to an underlying geometry and no contin- 
uous limit. Criticality is what links combinatorics and geometry, and 
what gives a meaning to approximation theorems. 

Definition 4.1. A discrete conformal map (<), p) is critical if there 
exists a discrete holomorphic map Z such that the quadrilateral faces 
<C>2 can be simultaneously embedded into rhombi in the complex plane. 

Because of the Gauss-Bonnet theorem, it is not possible to globally 
embed a compact surface into the plane, therefore we allow for an atlas 
of local critical maps with a finite number of fixed local conic singu- 
larities. When a continuous limit is taken, their number, angle and 
position should not change, and the theorem of isolated singularities 
helps us wipe them out as inessential. 

It is a simple calculation to check that if Z is a critical map, any 
discrete holomorphic function / G fi(<0>) gi yes r i se ; through (12.51) to a 
holomorphic 1-form fdZ. 

For a holomorphic function /, the equality / dZ = is equivalent to 
/ = Xe for some A G C with e the biconstant e{T) = +1, e(T*) = —1. 

Following Duffin [TH [12], we introduce the 

Definition 4.2. For a holomorphic function f, define on aflat simply 
connected map U the holomorphic functions f', the dual of f, and f, 
the derivative of f , by the following formulae: 

(4-1) /t(*) :=*(*)/(*), 

where f denotes the complex conjugate, e = ±1 is the biconstant, and 

(4.2) f'(z):=^(£pdZ^ +\e, 

defined up to e. 

It is an immediate calculation [I] to check the following 
Proposition 4.3. The derivative f fulfills 

(4.3) df = fdZ. 

4.2. d operator. For holomorphic or anti-holomorphic functions, df 
is, locally on each pair of dual diagonals, proportional to dZ, resp. dZ, 
we define d and d operators (not to be confused with the boundary 
operator on chains) that decompose the exterior derivative into holo- 
morphic and anti-holomorphic parts yielding 

dfAdf= (\df\ 2 + \df\ 2 )dZAdZ 
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where the derivatives naturally live on faces: 

In the continuous theory, for any derivable function / on the complex 
plane, the derivatives d = 4- + i 4- and 5 = A with respect to 

r ' ax ay ax ay r 

z = x + iy and z = x — iy yield 

f(z + z ) = f(z ) + z(df)(z ) + z(5f)(z ) + o(\z\). 

These derivatives can be seen as a limit of a contour integral over a 
small loop 7 around z : 

(df)(zo) = lim — % — <b fdz, (df)(z ) = - lim — % — & fdZ, 

7-*o 2A{^) J 7 7^0 2.4(7) J-y 

which leads to the following definitions in the discrete setup: 

d:C°(0) C 2 (0) 

f » df= [{x,y,x',y') >-> 2A[x ; yiX ,, y , } j fdZ], 

d:C°{0) -> C 2 (0) 

f » df= [{x,y,x',y')^- 2Aix l xf>yf) j fdZ]. 

A holomorphic function / verifies df = and (with Z(w) noted simply 
u) 

f(y') - f(y) f(x') - f(x) 



df(x,y,x,y) 



y' — y x' — x 



Notice that the statement / = (J df dz) has no meaning, d is not a 
derivation endomorphism in the space of functions on the vertices of 
the double. 

On the other hand, these differential operators can be extended 
(see [6]) into operators (the Kasteleyn operator) d 20 , B 20 : C 2 (<0>) — > 
C°(<0>) simply by transposition, d 2 o = — *<9o2, leading to endomorphisms 
of C°(<0) © C 2 (<£>). They are such that their composition, restricted to 
the vertices <0>o, gives back the laplacian: 

A = ±(dod + 5od). 

Furthermore, the double derivative d 20 o d 02 is a well defined endo- 
morphism of C°(<0>). 
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4.3. Discrete exponential. 

Definition 4.4. For a constant A G C, the discrete exponential exp(:A: Z) 
is the solution of 

exp(:A:0) = 1 

(4.4) dexp(:X:Z) = A exp(:A: Z) dZ. 

We define its derivatives with respect to the continuous parameter A: 

(4.5) Z- k - exp(:X: Z) := _ e xp(:A:Z). 

The discrete exponential on the square lattice was defined by Lelong- 
Ferrand [13] , generalized in [TJ] and studied independently in P, [15] . 
For |A| 7^ 2/5, an immediate check shows that it is a rational fraction 
in A at every point: For the vertex x = $ e > 

(4.6) -p(:A:x) = n rT i-^ 

k i 

where (6k) are the angles defining (5e l9k ), the set of (Z-images of) <0>- 
edges between x and the origin. Because the map is critical, Eq. (14.61) 
only depends on the end points (0,x). It is a generalization of a well 
known formula, in a slightly better version, 

, . , r Xx\ n ^X 2 x\ fl + tY 

4.7 expAx = 1 + — +0 = ^ + O — ^- 

\ n J n \^~^J 

to the case when the path from the origin to the point x = Yli ^ = 
^2Se ldk is not restricted to straight equal segments but to a general 
path of 0(\x\/5) segments of any directions. 

The integration with respect to A gives an interesting analogue of 
Z'~ k - exp(:A:Z). It is defined up to a globally defined discrete holo- 
morphic map. One way to fix it is to integrate from a given Ao of 
modulus 2/5, which is not a pole of the rational fraction, along a path 
that doesn't cross the circle of radius 2/5 again. 

Proposition 4.5. For point-wise multiplication, at every point x G (}q, 
(4.8) exp(:A: x) ■ exp(: — A: x) — 1. 

The specialization at X = defines monomials: 



(4.9) 



Z- k - = Z :fc: exp(:A:Z)| A=0 
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which fulfill Z :k - = k f Z'- k 1: dZ . The anti-linear duality t maps ex- 
ponentials to exponentials: 

(4.10) exp(:A:) t = exp(:^:). 

In particular, exp(:oo:) = 1' = £ is the biconstant. 

Proof \4~l>[ The first assertion is immediate. 

The derivation of (14.41) with respect to A yields 
(4.11) 

d k ( d k d k ~ l \ 

d g^j- exp(:A: Z) = I A exp(:A: Z) + k -^-^ exp(:A: Z) j dZ 

which implies ( 14.9ft . 
Derivation of exp(:A:)t gives, 

(4.12) (exp(:A:)t) / = 1 (£ exp(:A:) dZ^j + fie 

4 /exp(:A:) - l\ f 

+ fie 



5 2 V A 
= ^ exp(:A:) t + ve 

with /i, i/ some constants, so that the initial condition exp(:A:0)^ = 1 
at the origin and the difference equation d exp(:A:)t = exp(:A:)^ dZ 
yields the result. □ 

Note that it is natural to define exp(:A: (x — Xq)) := ^r^Q as a 
function of x with xq a fixed vertex. It is simply a change of origin. 
But apart on a lattice where addition of vertices or multiplication by 
an integer can be given a meaning as maps of the lattice, there is no 
easy way to generalize this construction to other discrete holomorphic 
functions such as exp(:A: (x + ny)) with x, y £ <C>o and n £ Z. 



4.4. Series. The series YlkLo A § ' > wherever it is absolutely conver- 
gent, coincide with the rational fraction ( 14.61) : Its value at the origin is 
1 and it fulfills the defining difference equation ( 14. 4ft . Using Eq. ( 14.91) . 
a Taylor expansion of exp(:A:x) at A = gives back the same result. 
We are now interested in the rate of growth of the monomials. 
Direct analysis gives an estimate of Z [k: : 
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Proposition 4.6. For x G <0>, at a combinatorial distance d(x, O) of 
the origin, and any k EN, 



(4.13) 



Z :k: (x) 



k\ 



< 



a 



a 



d{x,0) 



5 



a 



for any a > 1 arbitrarily close to 1. 



Corollary 4.7. The series Yl 



|A|<|- 



k=0 k\ 



is absolutely convergent for 



Proof U^j It is proved by double induction, on the degree k and on the 
combinatorial distance to the origin. 

For k = 0, it is valid for any x since = 1 + > 1, with equality 
only at the origin. 

Consider x G ^> a neighbor of the origin, = 5e j61 , then an 

immediate induction gives for k > 1, 



(4.14) 



fc! 



a+l 

Q-l 



a fc > 2. 



which fulfills the condition Eq. (I4.13P for any k > 1 because 
This was done merely for illustration purposes since it is sufficient to 
check that the condition holds at the origin, which it obviously does. 

Suppose the condition is satisfied for a vertex x up to degree k, and 
for its neighbor y, one edge further from the origin, up to degree k — 1. 
Then, 



(4.15) 



:fc:. 



■I' 



r-.k-V.I 



k\ k\ 

in absolute value fulfills 



(x) + Z-- k - 1 -jy)Z(y)-Z(x) 
(k-l)\ 2 



Z-- k \y) 



k\ 
(4.16) 



< 



a 



d(x,0) 



a — 1 

a + : 

a — 1 

a + 1 

a — 1 



fe-i 



a 



a 



(x.O) 



d(y,0) 



a 



thus proving the condition for y at degree k. It follows by induction 
that the condition holds at any point and any degree. □ 
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4.5. Basis. The discrete exponentials form a basis of discrete holo- 
morphic functions on a finite critical map: given any set of pair-wise 
different reals {A^} of the right dimension, the associated discrete expo- 
nentials will form a basis of the space of discrete holomorphic functions. 
See [8] for the formula 

(4.17) f(x) = -^Jg(X)exp(:X:x)dX 

for a certain fixed contour 7 in the space of parameters A, and the 
definition of g(X) as a fixed contour integral in <0> involving /. 

The polynomials however don't form a basis in general: the combi- 
natorial surface has to fulfill a certain condition called "combinatorial 
convexity" . A quadrilateral, when traversed from one side to its oppo- 
site, define a unique chain of quadrilaterals, that we call a "train-track" . 
The condition we ask is that two different train-tracks have different 
slopes. 

On a combinatorially convex set, the discrete polynomials form a 
basis as well. 

4.6. Continuous limit. In a critical map, where quadrilaterals are 
mapped to rhombi of side S, identifying a vertex x with its image Z(x). 

The combinatorial distance d$(x,0) is related to the modulus \x\ 
through 

(4.18) ^(x,0)^<M<^ (X)0 ) 

where 8 m is the minimum of all rhombi angles. When the rhombi 
don't flatten, the combinatorial distance and the modulus (over 5) are 
equivalent distances. 

Lemma 4.8. Let (ABCD) be a four sided polygon of the Euclidean 
plane such that its diagonals are orthogonal and the vertices angles are 
in [77, 27T — 77] with 77 > 0. Let (M, M') be a pair of points on the polygon. 
There exists a path on (ABCD) from M to M' of minimal length I. 
Then 

MM' sin 77 

i ~ ~^T' 

It is a straightforward study of a several variables function. If the 
two points are on the same side, MM' = t and sin 77 < 1. If they are 
on adjacent sides, the extremal position with MM' fixed is when the 
triangle MM'P, with P the vertex of (ABCD) between them, is isocel. 
The angle in P being less than 77, MAL > s i n 2 > ££2Z_ jf the points 
are on opposite sides, the extremal configuration is given by Fig. 02., 
where = ^.D 
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1. M, M' on adjacent sides. 2. M, M' on opposite sides. 

Figure 8. The two extremal positions. 

A function / : ^> ~~ ¥ C on the combinatorial surface can be extended 
to a function on the image of the critical map in the complex plane 
/ : U — ► C by stating that / (Z(x)) = f(x) for the image of a vertex, 
and extend it linearly on the segments [Z(x), Z(y)} image of an edge, 
and harmonically inside each rhombus. 

Theorem 4.1. Let (<0>fc) a sequence of simply connected critical maps, 
U the non empty intersection of their images in the complex plane and 
a holomorphic function f : U — > C. If the sequence of minimum angles 
are bounded away from and the sequence of rhombi side lengths (5k) 
converge to 0, then the function f can be approximated by a sequence 
of discrete holomorphic functions f n G fl(()k) converging to f . The 
convergence is not only pointwise but C°° on the intersection of im- 
ages. Conversely a converging sequence of discrete holomorphic func- 
tions converges to a continuous holomorphic function, in particular the 
discrete polynomials and discrete exponentials with fixed parameters. 

Corollary 4.9. On a Riemann surface, any 1-form can be approx- 
imated by a sequence of discrete holomorphic 1- forms on a refining 
sequence of critical maps with fixed conic singularities. 

The proof relies on the convergence of polynomials seen as iterated 
primitives of the constant function. 

Lemma 4.10. Given a sequence of discrete holomorphic functions (fk) 
on a refining sequence of critical maps, converging to a holomorphic 
function f, the sequence of primitives (J fkdZ) converges to f f(z) dz. 
Moreover, in the compact case, if the convergence of the functions is of 
order 0(51), & s ^ a V s this way for the primitives. 

Proof \4M% Suppose that we are given a sequence of flat vertices 
Ok G ()k where the face containing the fixed flat origin O G U is 
adjacent to Ok- For a given integer k, let F k the continuous piecewise 
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harmonic extension of the discrete primitive f Q fk dZ to U. We want 
to prove that for any x G U, the following sequence tends to zero 



(4.19) 



(F k (x) - F k {0)) 



f(z) dz 



o 



fceN 



For each integer k consider a vertex Xk G <0>o on the boundary of the 
face of ()2 containing x. 

We decompose the difference ( 14.191) into three parts, inside the face 
containing the origin O and its neighbor O^, similarly for x and Xk, 
and purely along the edges of the graph <0 fc itself. 



F k (0)) 


- r nz)dz 






Jo 




(F k (x) 


- F k (xk)) + 


< 


F k {x) - 


-Fk(xk)- [ 

Jx 


(4.20 


fk(O k ) 


-h{o)- I 

J 



•o* 



f(z) dz 



o 



+ 



J-k 



fkdZ 



eh 



f(z) dz 



f(z) dz 



o 



On the face of <0> containing x, the primitive £ i— > J x f(z) dz is a holo- 



morphic, hence harmonic function as well as £ i— > Fk(C,). By the maxi- 
mum principle, the harmonic function £ i— > -Ffc(£) —Fk(xk) — f5 /(z) <iz 
reaches its maximum on that face, along its boundary. The difference 
of the discrete primitive along the boundary edge (x k ,y) G Oi at the 
point £ = (1 — A)xfe + Ay is equal by definition to 

(4.21) F fc ((l - \)x k + Xy)- F k (x k ) = X(y - x k ) fkM + h{y) . 

The holomorphic / is differentiable with a bounded derivative on [/, 
so averaging the first order expansions at Xk and y, we get 



f(z)dz = X(y-x k ) 



X(y - x k ) 



f(x k )+f(y) + {y _ ^ )2 A 2 f(x fc ) + (l-A) 2 f(y) 
+o((£-x fc ) 3 )+o((£-y) 3 ) 



(4.22) 
therefore 

(4.23) 



|.F fc (x) - F k (x k ) 



f(z)dz\ = 0(S 2 k ). 



■n,: 
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Similarly for the term around the origin. 

By definition of fk, the 1-form fk(z) dz along edges of the graph <0> 
is equal to the discrete form fkdZ so that f@ k fk dZ = fk(z) dz on 
a path along (} edges. Therefore the difference 



(4.24) 



x k r x k 

f k dZ- / f(z)dz 

O k JOk 



< 



o k 



f k {z) - f{z)) dz 



is of the same order as the difference \fk{z) — f(z)\ times the length 
£(jk) °f a path on <0>fc from O fc to Xf.. This length is bounded as i^k) < 
. 4 e \xk — Ok\- Since we are interested in the compact case, this length 
is bounded uniformly and the difference (14.241) is of the same order as 
the point- wise difference. We conclude that the sequence of discrete 
primitives converges to the continuous primitive and if the limit for 
the functions was of order 0(5 2 ), it remains of that order. (} 

The discrete polynomials of degree less than three agree point-wise 
with their continuous counterpart, Z :2: (x) = Z(x) 2 . 

A simple induction then gives the following 

Corollary 4.11. The discrete polynomials converge to the continuous 
ones, the limit is of order 0(5|). 

Which implies the main theorem: 

Proof \4^1\ On the simply connected compact set U, a holomorphic 
function / can be written, in a local map z as a series, 

(4.25) = 

fceN 

Therefore, by a diagonal procedure, there exists an increasing integer 
sequence (N(n)) nGN such that the sequence of discrete holomorphic 
polynomials converge to the continuous series and the convergence is 
C°°. 



(4.26) > ^ a k Z- h -- | - /. 








4.7. Cross-ratio preserving maps. Once the isometry Z is cho- 
sen, holomorphicity of a function / can be written on a quadrilateral 

(x, y, x', y') G <^>2, writing x = Z{x) for a vertex x G <)o, as 

f(y') - f(y) y' -y 



(4.27) 



f(x') — f(x) x' — x 
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and / is understood as a diagonal ratio preserving map, and each value 
at a corner vertex can be linearly solved in terms of the three others. 

A quadratic version is given by the cross-ratio preserving maps: A 
function / is said to be quadratic holomorphic iff 

U9H] u\y) - m) uw) - m) _ (?-*) w - ^ 

{ ' ) (/(*) - f(y')) (f(x') - f(y)) (x - y>) {x> - y) ■ 

A rhombic tiling gives rise to two sets of isoradial circle patterns: a 
set of circles of common radius 5, whose centers are the vertices of T 
and intersections are the vertices of T* and vice-versa. Two interesting 
families of cross-ratio preserving maps are given by circle patterns with 
the same combinatorics and intersection angles as one of these two circle 
patterns. 




Figure 9. A circle pattern with prescribed angles as a 
cross ratio preserving map. 




Figure 10. The cross-ratio q is given by the intersection angles. 

A change of coordinates helps understanding diagonal ratio preserv- 
ing maps as a linearized version of the cross-ratio preserving maps. We 
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will say that the function w : ^>2 — ► C solves the Hirota system if, 

around a face (x, y, x', y') G <^2> 

(4.29) 

(y— x)w(x)w(y)+(x— y')w(y')w(x)+(y'— x')w(x')w(y')+(x' — y)w(y)w(x') = 0. 

This is to be understood as a quadratic version of the Morera theorem 
§ f dz = and w is a half of the derivative of a holomorphic function: 

Proposition 4.12. If w solves the Hirota system, then the function 
f : <J)2 —> C defined up to an additive constant by 

(4.30) /(y) - f(x) = {y- x)w{x)w{y) 

is quadratic holomorphic. 

Proof \4^T^ The function / is well defined because the associated 
1-form is closed by definition of the Hirota system. The function w 
disappears in the cross-ratio of /, leaving the original cross-ratio. <£> 
Conversely, a quadratic holomorphic function defines a solution to the 
Hirota system, unique up to multiplication by A on T, 1/A on T*. Con- 
cerning circle patterns families, w is real on the centers and unitary on 
the intersections, and encodes the variation of radius, resp. of direction 
of the image of the circle: 

(4.31) f(y)-f(x)=r(x)e ie ^( V -x). 

Proposition 4.13. The logarithmic derivative of the Hirota system 
associated to a family of cross-ratio preserving maps is a diagonal ratio 
preserving map. 

In other words, for (1 + eg)w to continue solving the Hirota system 
at first order, the deformation g must satisfy 

( 4 32 ) g(y') - gjy) = fit/) - f{y) 

g(x')-g(x) f(x')-f(x)' 
Proof 4-13, The e contribution of the closeness condition (I4.29p for 



(1 + eg)w gives 

(g(x) + g(y)) f(x) f(y) (y - x) + {g{y) + g{x')) f(y) f(x') (x' - y) 
+ (g(x) + g(y')) w(x') w(y') (y' - x') + (g(y) + g(x)) w{y') w(x) (x 

which reads, refering to /: 



g(x') - g(x) f(x') - f(x) 
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4.8. Baecklund transformation. The Baecklund transformation is 
a way to associate, to a given solution of an integrable problem, a family 
of deformed solutions. The two problems under consideration here are 
the linear and quadratic holomorphicity constraints on each face. They 
are given by a linear, resp. quadratic algebraic relation between the 
four values of a solution at the vertices of each face. These relations 
involve only values supported by the edges of the rhombus, which are 
equal on opposite sides, namely the complex label y — x — x' — y'. 

The Baecklund transformation is defined by imposing such constraints 
over new virtual faces added over each edge, with "vertical edges" la- 
belled by a complex constant A: 

Definition 4.14. Given a linear holomorphic function f G ^(<)), com- 
plex numbers u,X G C, its Baecklund transformation f\ = B%(f) is 
defined by 



(4.33) 



/a(0) = u, 
f x (x) - f(y) X + x 



fx(y) ~ f{x) X + y-x 



Given a quadratic holomorphic function f , complex numbers u, X G C, 
its Baecklund transformation f\ = B%(f) is defined by 



(4.34) 



/a(0) = u, 
fx(y)-fx(x) f(x)-f(y) (y - xf 
h{x) - f(x) f{y) - f x {y) A 2 



The right hand sides are the values respectively of the diagonal ratio 
and cross-ratio of a parallelogram faces of sides (y — x) and A seen as 
"over" the edge (x, y) G 




Figure 11. The face (x x , yx, x' x , y' x ) "over" the face 
(x,y,x',y') G <^ 2 - 



DISCRETE RIEMANN SURFACES 



33 



Proposition 4.15. This transformation is well defined in the critical 
case. 

This condition, called three dimensional consistency is an overdeter- 
mination constraint: if the cube "over" the face (x,y,x',y') G {>2 is 
split into two hexagons along the cycle (y, x', y', y' x , X\, y\), one can see 
that, given values at these six vertices, the values at the centers of each 
hexagons, namely at x' x and at x are overdetermined. 




Figure 12. The cube split into two hexagons yielding 
equivalent compatibility constraints. 



Therefore only certain values at the six vertices are allowed, defined 
by two algebraic relations between them. The compatibility condition 
is that these two algebraic relations are equivalent. It is a simple com- 
putation to check it is the case for critical maps. 

This transformation verifies 

(4.35) B'V> (Bt(f)) = f 

for any (u, A). It is an analytic transformation in all the parameters 
therefore its derivative is a linear map between the tangent spaces, that 
is to say between diagonal ratio preserving maps, 

(4.36) dB u x (f):Q(f) -> (£*(/)). 

It is not injective and I define the discrete exponential at / as being 
the direction of this 1-dimensional kernel. It can be characterized as a 
derivative with respect to the initial value at the origin: 

(4.37) exp u (:A:/) := | u=/(0) G ker (dB%(f)) 

because B x (B x _x(g)) = g for all A, g and v. 

As in the discrete exponential case, the value of the Baecklund trans- 
formation at a given vertex is the image of values at neighbouring 
vertices by a homography. These homographies can be encoded as 
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9 = B u x (f) 



Figure 13. The discrete exponential exp u (:A:/) is the 
kernel of the linear transformation dB'^(f) (here u = 
f(0))- 



projective operators L(e; A) G GL 2 (C)[X] lying on the edges e G (} 2 , 

called a zero curvature representation: 

(4.38) 

'X + y-x -2(y-x)(f(x) + f(y)f 
A + x — y 



L((x,y); A) 



for the linear case, 



(4.39) 
L((x,y); A) 



1 -(y-x)w(y) 
-X(y — x)/w(x) w(y)/w(x) 



for the Hirota system. 



Then we define [16] the moving frame \P : <C>2 — > GLi2(C)(\) by a pre- 
scribed value at the origin and recursively by $(y; A) = L((x, y); X)^(x; A) 
and its logarithmic derivative with respect to A 

d*(e;A), 



(4.40) 



A{e-X) 



dX 



is meromorphic in A for each edge e. We call /, resp. w isomonodromic 
if the positions and orders of the poles don't depend on the edge e. The 
two points discrete Green function (the discrete logarithm) G(0,x), 
inverse of the Laplacian in the sense that 



(4.41) 



AG{Q, 



5, 



can be constructed as the unique isomonodromic solution with some 
prescribed data [8], which allows us to give an explicit formula for it, 
recovering results of Kenyon [15] : an integral over a loop in the space 
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of discrete exponentials: 

(4.42) G(0, x) = -— ^ / exp(:A: x) !^-^ dX 

O 7T % Jq A 

where the integration contour C contains all the possible poles of the 
rational fraction exp(:A:x) but avoids the half line through —%. It is 
real (negative) on half of the vertices and imaginary on the others. 
Because of the logarithm, this imaginary part is multivalued. 
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